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CN ■ Abstract 

■ We explicitly show that the one loop IR correction to the two-point function in de Sitter 

■ space scalar QFT does not reduce just to the mass renormahzation. The proper interpretation of 

■ the loop corrections is via particle creation revealing itself through the generation of the quantum 
averages (a^ Cp), (apa_p) and {apU'tp), which slowly change in time. We show that this observation 
in particular means that loop corrections to correlation functions in de Sitter space can not be 

pi I obtained via analytical continuation of those calculated on the sphere. 

^ I We find harmonics for which the particle number (a^ Op) dominates over the anomalous expec- 

tation values {ttpQ-p) and {a^a^p). For these harmonics the Dyson-Schwinger equation reduces 
\ in the IR limit to the kinetic equation. We solve the latter equation, which allows us to sum up 

\ all loop leading IR contributions to the Whiteman function. We perform the calculation for the 



in 
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principle series real scalar fields both in expanding and contracting Poincare patches. 



INTRODUCTION 



There are large infrared (IR) loop contributions to the correlation functions even due to massive 
^ ! fields in de Sitter (dS) space. The question is whether after the summation over all leading 

5^ , IR loop contributions the correlation functions will be finite or will grow unboundedly. If the 

correlation functions do grow in the IR limit, then this is the sign that back-reaction of the 
quantum fluctuations on the background dS geometry is not negligible. The literature on this 
subject is vast, see e.g. []]|-j2o|. 

In general to sum the loop contribution one has to solve the Dyson-Schwinger equation (DSE). 
However, the interest is in the solution of this equation in the extreme IR limit, i.e. in the sum 
of the leading IR corrections. Colloquially speaking in the latter limit quantum coherence is lost 
and the legs in the loops sit on mass-shell. The situation is analogous to the one discussed in 



a similar context e.g. in [21[-23]. As the result, known in condensed matter physics, the DSE 



for non-stationary diagrammatic technic reduces to the Boltzmann equation (see e.g. [2j], [251]). 
In other words the classical Boltzmann kinetic equation (KE) allows to sum up the leading IR 
corrections. We concisely review the relation between DSE and KE in the Appendix. 

In this note we derive the Boltzmann equation in the Poncare patch (PP) of dS space. For the 
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earlier papers on the KE in curved spaces in general and in dS space in particular see e.g. |26l|.[27|. 
Unlike the equations of those papers our equation does not really describe the dynamics of the 
occupation numbers of particles in dS space, because in its derivation we use the exact harmonics on 
dS background rather than the plain waves. It describes the dynamics of the occupation numbers 
of the waves with wavelength comparable to the horizon size. We would like to study the physics 
due to long wavelength fluctuations because IR effects are dominated by such objects. Still we call 
them as particles throughout the paper. 

Our eventual goal is to see whether the back-reaction from the quantum fluctuations in dS space 
is negligible or not. But for the beginning we neglect back-reaction and consider self-interacting 
massive QFT on the fixed dS background. We restrict our attention to the behavior of the solution 
of the KE in the extreme IR limit. 

The consideration of the exact harmonics in the background fields brings several complications. 
The main one is generic for the curved space-times and is due to that even for the fundamental 
questions in curved space-times it is necessary to specify boundary conditions, so to say. This is 
because there is no any preferable reference frame in a general curved space-time. 

In dS space this complication reveals itself through the crucial difference between Cauchy prob- 
lems for the KE as defined in global dS and in PP. The main difference is due to the different 
geometry of the Cauchy surfaces in these two situations. The PP of dS space is the grey region in 
the rectangle shown on the figure 1 — Penrose diagram of the dS space. The left and right sides 
of the rectangle are glued to each other. The Cauchy surfaces in the PP are depicted as the lines 
close to the hyperbolas shown in this figure. As we go to the past infinity the Cauchy surfaces 
degenerate to the boundary of the PP. The Bunch-Davies (BD) vacuum means that there are no 
positive energy states on the boundary of the PP. Complementary to the grey part of the diagram 
is the contracting PP. 
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The Cauchy surfaces in global dS space are depicted by the lines of the type shown on the figure 
2. E.g. the Euclidian vacuum means that there are no positive energy states on such a surface 
going nearby the neck of dS space. In such a situation we have a different state at the boundary 
of the PP. Thus, even despite that the Whiteman function for BD vacuum state in PP coincides 
with the one for Euclidian vacuum in global dS, the QFT dynamics in both situations can be quite 
different. In particular, global dS space contains as well the contracting PP. We point out that 
the behavior of solutions of the KE in the contracting patch is quite different from the expanding 
one. Furthermore, from the cosmological point of view one needs to study global dS because it sets 
up initial conditions for the PP — for the inflation. But this issue demands a separate study and 
will not be addressed in the present paper. Here we restrict our attention to the expanding and 
contracting PP of dS space. 

Another complication is peculiar for homogeneous time-dependent backgrounds with horizons 
such as dS space. It comes from the fact that unlike Minkowski and AdS space there is no unique 
choice of harmonic basis which diagonalizes free Hamiltonian once and forever. To address this 
complication we do not specify harmonics in the general formulas. We attempt to use formulas 
which have a generic application in dS space or even for more general FRW space-times with flat 
spacial sections. At the end we of cause study physics due to different specific choices of harmonics. 

Yet another complication is as follows. The regular KE for particles describes the change in time 
of the particle density Up = {at a^) . At the same time one assumes that the anomalous quantum 
averages Kp = (ap-O-p-) and k* = {at a'^t) are vanishingly small. In background fields in general 
this is not the case any more. Only if one finds the appropriate equilibrium state he can neglect 
Kp. In the paper the KE of the type interesting for us was derived. However in that paper the 
analysis of the presence of the anomalous quantum averages Kp and k* was not done. As the result 
the wrong choice of the harmonics was made there. 
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Furthermore, waggly speaking KE is valid if its solution Hp is slow function of time in comparison 
with the corresponding harmonics. The consideration of the long wavelength fluctuations poses 
the question about the existence of such a separation of scales in KE. However, we check whether 
Hp is slow enough or not for every explicit solution of the equation, which we find. 

Finally, in this paper we discuss the scalar fields from the principle series, which have masses 
m > 3/2 in units of the dS curvature. The dynamics of the light fields, m < 3/2, i.e. from the 
complementary series, can be quite different. It goes without saying about the massless fields. 

II. SPECIFICATION OF HARMONICS 

The goal of this paper is to derive kinetic equation (KE) for the exact harmonics in de Sitter 
(dS) space. To do that we have to specify which harmonics one should use. As we explain in the 
next section the proper KE can be derived only in such circumstances when one can neglect (a^ a_p) 
in comparison with (0^0^), where a and a~*~ are annihilation and creation operators. Below we 
will be able to specify such harmonics for which {ap fl-p) is negligible in comparison with (a^ ap) 
at the low physical momenta prj — t- 0, where rj is the conformal time in the dS spacc^. 

The way we will do that is as follows. The Kcldysh propagator carries statistical information 
about the theory and the state in question (see the discussion below). The calculation of the one 
loop contribution to the Keldysh propagator will allow us to estimate the behavior of (a+Op) and 
(opO-p). For some harmonics {Up ap) will be of the same order as (apa-p), while for the others 
[apQ-p) will be suppressed in the IR limit. 

We would like to consider scalar field theory in the background space-time with the metric: 

ds'^ = a^{r]) [dr]^ -dx^] . (1) 

In this paper we always consider (+,—,—,—) signature. Although all our formulas can be straight- 
forwardly generalized to arbitrary dimensions we restrict ourselves to 4D. 

The choice 0(77) = 1/r], < rj = e~* < +00 corresponds to the PP of dS space. We put the 
Hubble constant to one H = 1. Past infinity of the PP corresponds to 77 = e~* — t- +00, while 
the future infinity is at 77 = e~* — t- 0. Wc prefer the definition t] = e^* > instead of the more 
standard one rj = — e^* < to avoid dealing with non-integer powers of the negative quantity. The 
contracting PP corresponds as well to 0(77) = 1/r], but now 77 = e* and past infinity corresponds to 
77 = 0, while future infinity — to ry = +00. So in contracting PP conformal time flows in the reverse 
direction with respect to (wrt) the expanding PP. Below we always talk about either expanding or 
contracting PP of dS space, but prefer to keep generic 0(77) in those formulas which may be applied 
to more general space-times. 

We consider the theory of the real massive scalar particle with the cubic self-interaction: 



^ The extreme IR limit for the physical momentum corresponds to the future infinity for the fixed comoving mo- 
mentum in the expanding Poincare patch (PP) of dS space. 



5 



S= j d^x j dr]^/\g\C[(t)\= j d^x j dr}a^{r]) 

J J OO J J OO 



2^ (5m'/') _ 



(2) 



Although the cubic potential has a run away instability we have chosen it to simplify all our 
formulas. A short comment on how the instability of the cubic theory reveals itself in the KE can 
be found in the first part of the Appendix. Conceptually all our calculations are valid as well for 
the theory with (/)^ self-interaction term. In fact, it is not difficult to write the answer for the latter 
theory once the answer for the cubic self-interaction is known. 

A. Tree— level two— point function 

To set up the notations wc have to start with the tree-level two-point function. The expansion 
of the free field in terms of the normalized harmonics is 



0(7?, x) = j 



d^k 



where (pkiVi = 9k{v) ^ ^"^^ some solutions of the Klein-Gordon equation. 



(3) 
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(p{v,x) = 0, 



(4) 



in the metric under consideration. Concretely gkiv) = V'^^'^ h{kri) / \/2, where h{krj) is some properly 
normalized solution of the Bcsscl equation, a(ry) = l/rj. 

From the energy momentum tensor T^,^ = d,j_(f)du4> ~ Qi^u ^[4>\j we find the Hamiltonian: H(r]) = 
(rj) J d^x Too (77) . The free Hamiltonian looks as^ : 

Hoiri) = J d^k [a'^ ^^(r/) + a_jt Bk{rj) + h.c] , 




+ [k' + a\v)m']gl\. 



(5) 



The harmonics, which simultaneously solve the Klein-Gordon equation and Bk = can be found 
only in those regimes, when the background field is switched off and, hence, and Bj^ are time 
independent. Then the harmonics just coincide with the plain-waves. In general the solutions of 
the equation B^ = and of the Klein-Gordon one do not coincide, because B^ is time dependent. 



^ After the normal ordering of the free Hamiltonian one has to absorb the infinite vacuum energy into the redefinition 
of the cosmological constant. Such a quantity is time independent if expressed in terms of the physical momentum 
p/a{ri) and attributed to the Hamiltonian conjugate to the time t rather than to the conformal time r) = e~*. 
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Hence, the free Hamiltonian is not diagonal. Furthermore, one can not choose the sohitions of the 
equation Bf^ = as g^, because then the corresponding and do not obey the Heisenberg 
algebra. 

We do not specify harmonics in the formulas which have general application, but below we are 
going to consider particular cases. Let us make a few comments about the standard choice of the 
Bunch-Davies (BD) harmonics in the PP 13]: 



^=^^2--, (6) 
where Til^' (x) is the Hankel function. We choose Til^' (x) as the positive energy harmonics instead 

(2) 

of Til^ (x) because of the reverse order of the time flow: oo — )• — )• 0. It is not hard to see that 
the harmonics in question diagonalize the Hamiltonian at the past infinity — )• oo, because in this 
limit gkif]) ~^ '^^^^^ ~^ ^' ^^^^^ ~^ k. In the future infinity — >• there is a solution 

of the Klein-Gordon equation, which only approximately solves = for the very massive fields 
m > 3/2. 

The Keldysh propagator is defined as'^ {rii,ri2,\x — y\) = ^ {{<p{rii,x), (f>{ri2,y)}) ■ Due to 
the spacial homogeneity of dS space we find it more convenient to make the Fourier transform 
along the spacial directions — -0^(7/1,172) = / cfir {7]i,r]2,r) e~^^^ . Furthermore, to simplify 
the formulas we define d^ijprji, pr]2) = (??i??2)^^ -0^(r/i, 772). At tree-level this quantity is equal 
to: 



9kiv) = 



do{pVi,PV2) = ^ [h{pvi) h*{pri2) + h*{prii) h{pr]2)] (l + 2 (a+ ap)) + 

+KpVi) Kvm) («p a-p) + h*{pr]i) h*{pri2) (a+ atp) , (7) 

if we average wrt an arbitrary state. If the average is taken wrt the vacuum state ap\) =0, then 
(apG^p) = dip) = 0, and Op) = 0. Different choices of the harmonics h{pr]) correspond to 
the different choices of the vacuum state. Using Gradshtein and Rizhik eq. 6.672.2, one can show 
that for the arbitrary choice of the solution of the Bessel equation h{pri) the Keldysh propagator 
is equal to: 



D''{m,V2,\x-y\) = Ci{z^-l) P\^^^{z) + C2{z'-1) ' QIi^^^{z), (8) 

where and Ql are associated Legendre functions; z = 1 + ^-q^ -q!^ is the hyperbolic 
distance between the two points in question on dS space; C1.2 are some complex constants whose 



^ Note that here we use a shghtly different definition of the diagrammatic rules than in the Appendix. The difference 
is in the Keldysh rotation — instead of the one used in the Appendix here we use i^^i = (<?!'+ + 4>-)/2, (j>q = 4'+ — ■ 
The rules in the main body of the text agree e.g. with [29|, while those used in the Appendix agree with [25l |. 
That is because here we would like to compare our formulas to [2^, while the formulas in the Appendix should be 
compared to [25I ]. 
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values depend on the particular choice of the Harmonics. E.g. for the BD harmonics ^ C2 = 
and the two-point function coincides with the one following from the analytical continuation from 
the sphere and having proper Hadamard behavior. 



B. One loop two point function 



In this subsection we basically repeat and generalize the calculations of [14} | and . We would 
like to find the leading one loop contribution to d^{pr]i, prj2) as and Tyi/r/2 = const 

when we start from the vacuum state at the past infinity. Using Schwinger-Keldysh diagrammatic 
technic one obtains the leading answer: 



df{pVi,Pm) ~ 7; KPm) h*{pri2) X 



1 



27r2 



2 7r2 



^/'i dk 



1/'' dk 
T 

I 



dxi dx2 {xi X2) 2 h 



p 
Ik 



dx2 (xi X2) 2 h* 



P 



Xl 



Xl 



h* 



p 



■X2 



h\xi) [h*{x2)Y- 
-h{pr]i) h{pr]2) x 



P 



X2 



h\xi) [h*{x2)f + C.C. 



(9) 



where under the k i nteg ral it is assumed that 1/?? = l/\/^?i^2 ^ k ^ p. This formula reduces to 



141 ] when iji = r]2- 



the one obtained in 

From the last expression we see that although we have started from the vacuum state, where 
(op a_p) = and (a+ Cp) = 0, these quantities are generated at one loop level. This can be traced 
back to the pair creation in dS space. Indeed in the course of evolution towards future infinity 
r/ — )• the density of the created particles appears to be: 



npir]) = (fflp ap) 



A2 
47r2 



1/'' dk 



'JT 1 1 dxi dx2 {xi X2)'^ h 



P 

-k^^y 



At the same time the anomalous quantum average is given by 



P 



h\xi) [h*{x2)r.m 



an d- 



2 7r2 



^/"^ dk 
T 



Xl 



dx2 {xi X2) 2 h* 



k 



Xl 



h* 



k 



X2 



h\xi) [h*{x2)]\n) 



As the side remark let us point out that if one by mistake were using the standard Feynman 
diagrammatic technic in the circumstances under consideration, he would never see the appearance 
of the terms proportional to h{pr]i) h{pr]2) or h*{pr]i) h*(p7]2) in loop contributions to the Whiteman 
function^. It is not hard to see that only terms proportional to h* [prji) h{pr]2) are generated in 
the loops within the standard Feynman technic. As well in the stationary flat space case the 
Schwinger-Keldysh technic gives the same answer as the Feynman one, which reveals itself through 
the cancelation of all terms contributing to Kp. 



* From the Whiteman function one can construct Keldysh and/or Feynman propagator. 
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Now let us calculate the leading one loop IR contributions for various choices of harmonics 
h{pr]). To start we choose BD harmonics For this choice, the xi and X2 integrals in ([9]) rapidly 
converge and are saturated around xi^2 ~ Z^, because in this case h{x) ~ e*^' as x — )• oo. Hence, 
because p/k <C 1 we can Taylor expand the h{px/k) functions around zero, using their leading IR 
behavior 



h{x) A+ x'^" + A. x-'f" , x^O, 



A, 



IT e~ 



2^^+2 r (1 + i/u) sinh(7r/i) 



A. 



vr e ""2" 



r(l - i^) sinh(-7r/u) 



(12) 



Furthermore, expanding /i(p??i,2) for small and then keeping only leading IR terms in the xi 
and X2 integrals, one can find the sum of the tree level and one loop leading IR contributions 



coth (vr^) 
2/i 



s''' + A+A*__ {pvf'^ X 



x<!l + ^lo (— 
27r2/i °^ \pr] 



dxidx2 (xiX2)2 h (xi) [h*{x2)] x 

+ c.c, 



e{x^-x2) ["^Y -e{x2-x^)(^^ 

X2j \X2 



(13) 



where 9{xi — x?) is the Heviside ^-function and s = r]i/r]2, f] = ■s/r]irj2- In deriving this result we 



have used that 



dxx^+'^'h'^ix] 



-2-K 



(14) 



Our result (jl3p reduces to the one obtained in [3| when s = 1. Thus, we see that contributions 
to rip = ap) and Hp = {ap a^p) are of the same order. As we explain in the next section this 
means that BD harmonics are not suitable to write the KE only for Up. 

Let us consider a different choice of harmonics — e.g. so called Jost functions at future infinity 



h{x) 



vr 



sinh(7r/i) 



Ji^ (3^) ■ 



They behave as h{x) ~ x*'^, when x — )• 0. If one substitutes these harmonics into 
take into account that these harmonics behave at past infinity as 



(15) 

I, he should 



h{x) 



TT 



4 sinh(7r fi) x 



fl—l X 



X — )■ 00. 



(16) 



Hence, the xi and X2 integrals in ([9|) have contributions around infinity due to the interference 
terms between e*^ and e~*^. They do not converge fast enough: they are saturated in the vicinity 
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of px/k ^ fj, rather than at x ~ /u. Hence, naively one can not Taylor expand h{px/k) around zero 
inside the xi^2 integrals. However, let us subtract from and then add to h'^{x) and under 
the x-\ 2 integrals the values of the interference terms, . ^f/ '"s : 

-^1^ ^ ' 4 sinn(7r /^} x 



vr e 



+ 



vr e 



4 sinh(7r fi) x 



4 sinh(7r fi) x 4 sinh(7r fi) x 

and of its complex conjugate are saturated around 



Then the xi^2 integrals of h'^{x 
X ~ ^ and one can Taylor expand h{px/k) around zero inside the corresponding expressions. At 
the same time the contributions from the additional integrals of 4 sinh(7r ^) x suppressed in the 
IR limit. Indeed, due to extra powers of q the momentum integrals, dq, are not divergent in the 
limit pi] ^ 0. 

Thus, the leading IR contribution to the two-point function in this case is as follows: 



1 + 





/■o 1 . 


{-) 


/ dxx'2^'^^ 


\pvJ 


J 00 



h'^ix) 



vr e 



4 sinh(7r /i) x 



(17) 



We see that for this choice of the harmonics the particle number density rip dominates over Kp 
in the extreme IR limit. As we explain in the next section this means that the harmonics under 
consideration are suitable for the derivation of the KE in dS space at future infinity. 

C. de Sitter vs. Sphere QFT 

Before going on with KE let us make a few comments about the possibility to formulate the dS 
QFT via analytical continuation from the sphere [l3| . The IR limit of the tree-level propagator 
for the BD harmonics is: 



do {pvi,pm) 



coth (vr/i) 
2il 



(18) 



Mass renormalization + A/i would lead to the one loop contribution of the following form: 

coth (vr/i) 



di {prii^Pm) 



2/i 



■s'^'iAn log(s) + A+A*_ {prif^ 2iAfi log(p??) +c.c.. 



(19) 



Note that one expects the correction A^u to be complex [29(]. The last expression can not reduce 
to the actual result (|13p for any choice of A/x. In the light of what we have been saying in the 
sections above it is conceptually misleading to interpret the one loop contribution to the Whiteman 
or Keldysh propagator in dS space as the mass renormalization, because its proper interpretation 
is in terms of particle creation — in terms of slow functions Up and Kp of the average conformal 
time 7] = ^T]\ T]2- 

At the same time the one loop Feynman diagrammatic technic calculation on the sphere can 
lead only to the mass renormalization: the authors of [3] expected that after the appropriate 
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subtraction of the UV divergences the remaining finite part would coincide with the one obtained 
by the direct calculation in dS space. However, now one can see that this can not happen. So it 
is not correct to define QFT on dS space via analytical continuation from the sphere, even in the 
situations when the tree-level dS correlation functions are indeed obtained via such a continuation. 

Analytical continuation from the sphere may describe the stationary situation in dS, if any, but 
not the approach to the stationarity. The stationary situation corresponds to /Cp = and Hp being 
independent of rj. In fact, in the stationary situation the two-point correlation function depends 
only on the time difference (on r]i/r]2 in terms of the conformal time) rather than on both of the 
times (r/i and 772) independently. Such a situation can be achieved in dS space only in the extreme 
IR limit. Prom the one loop calculation above one may see that in BD state stationary situation can 
not happen in principle. However, for the Jost functions at future infinity the stationary situation 
may happen if the result of the summation of all loops — np{r]) — becomes somehow independent 
of time in the future infinity. 



III. THE PHYSICAL MEANING OF k, 



The one loop calculation of the previous section shows that the two-point function behaves as 



d {pr]i,Pil2) = 2 [HpVi) h*{pr]2) + h*{pr]i) h{pr]2)] [1 + 2np {^/r]rfn)] + 
+h{pr]i) h{pr]2) Kp (y/vm) + h*{pr]i) h*{pr]2) K*p {^/mm) , 



(20) 



when pr)i^2 0. Here np{r]) and Kp{rj) are slow functions of their argument in comparison with the 
harmonics h{pr]). Note that they are functions of time only, i.e. homogeneous in space. That is a 
natural situation in such a homogeneous space as dS, especially for the small values of the physical 
momenta. 

Such a situation allows to write a system of KE for both Up = y a^) and Kp = y {ap a^p). 
Here V = (5(0) is the comoving volume of the spatial sections. From the normalization of the 
harmonics it follows that then Up is the density per physical volume — per a^(r/) V — because Up 
appears in the expression for the total number density as follows: N = J d^pa~^{r])np{r]). Here 
p is the comoving momentum, while p/a{r]) is the physical one. The particle number density per 
comoving volume is defined as rip = np/a^ij]). 

The interaction Hamiltonian has the form (up to zero mode treatment, which is concisely 
discussed in the appendix): 



X |3(5 (-fei + ^2 + 



^3) 



Hint = ^ j d^ki d^k2 d% a\v) X 
all 9k2 9k3 iv) flfei ak2 ak3 + 9ki gl^ gls iv) aki a^^ 



+ 5(ki + k2 + k3^ gki9k2 9k3{rj)ak,ak^aks+gl,gl^9l^iv)aki(^k2^k3 }' (^1) 

Using this Hamiltonian one can find the evolution in time of {a~^a) and (aa) and proceed along 
the same lines as in the first part of the Appendix. But now we have to take into account, when 
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perform Wick contractions, that not only (^a^ ap/^ = np6{p — p'), but as well (^apap/) = Kp5{p + p') 
and (^a+ ap^ = Kp5{p + p'). 

With the use of the following matrixes: 



iv,(,,. m) = ( -.fUiM \ P = ( » M (22) 

V i^piVi) 9p{V2) np{r]i) gp{r]2) J \^ ^ J 

the resulting system of KE can be written in a compact form. The real equation has the form: 



dnp{r]) 
drj 



{rirf) 



Re Tr 



5[p-ki~k2) Cpk^kiiv) 



1 + N;) Nk, N^^ - N; (1 + AT, J (1 + N, 



k2) 



+26(h-k2-p) Ck.kM m^il + Nk,){l + Np) - {l + N*,jNk,Np (??',??') + 



+5 (p + ki + k2j Dpk.kM [{I + Np) (1 + Nk,) (1 + Nk^) - Np Nk, Nk, 
while the complex one is as follows: 



(^',^0 }, (23) 



dKp{r]) 



+2(5 (ki 



dv " ' ' J (2^)^ Jno ivv' 

p-ki- Cpk.kM [(1 + Np) (1 + Nk,) (1 + iVfe) - Np Nk, Nk, iv', V')+ 



+5 [p + ki + k2] D*pk^k2iv) 



[l + Np)N*k^N,,^ 



Np {1 + Nl) (l + iV4) 



iv ,v) + 

], (24) 



r/o is the moment of time when we switch on the interactions. The notation [N ^ P N] means that 
we have to add to the explicitly written expression the same one where every N is substituted by 
the product P N; Re Tr means that one has to take the real part and the trace of the expression 
following after these signs; at the same time {p — t- —p) Tr means that one has to take the trace and 
add to the expression following after these signs the same one with the exchange p — t- —p. Finally 

Ckikika = g*k^ gk2 9k-i and Dk^kik^, = dki gki dk^- 

To convert (j23p . (|24p into a tractable form one has to assume that np(rj') and Kp(rj') are slow 
functions in comparison with the harmonics. The situation is similar to the one discussed in the 
Appendix. Then the first argument of the matrix Nk(ri',ri') can be taken to be instead of rj' . 
Furthermore, explicit check shows that one can extend the limits of the 77' integration inside the 
collision integral (CI) to 77 —t- 0, ryo ~^ 00. This does not make the rj' and k integrals divergent. 

Expanding the terms under the CI in (|23p one will encounter the standard expressions which 
appear in the CI written in the Appendix: 
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[(1 + Up) rifei nfca - np (1 + n^J (1 + n^J] 
[rifc^ (1 + nfcj (1 +np) - (1 +nA:i)nfc, 



2 "^pJ ' 

[(1 +np) (1 (1 + 71^2) -npTifc^nfcj] (25) 

These three contributions have the following physical meaning. The first term describes the com- 
petition between two processes. One that the wave with the momentum p decays into two waves 
= P- This process, corresponding to the term rip (l + rifc^) (l+n^j), appears with the minus 
sign in the CI because it describes the loss of the wave with the momentum p. The inverse gain 
process, corresponding to the term (1 + np)nk^ n^^ with the plus sign, is that when two waves 
merge to create the wave with the momentum p. 

The second term in (j25p describes as well two competing processes. The first process is that 
the wave with the momentum p can merge together with another wave (with the momentum 
^2) to create a third one (with the momentum ki). This is the loss process. The inverse gain 
process happens when a wave (with the momentum ki) decays in to two, one of which is with the 
momentum p. The coefficient 2 in front of this term is just the combinatoric factor because we can 
exchange ki and k2. 

The third term as well describes two processes. The gain process is when three waves, one 
of which is with the momentum p, are created out of vacuum. The loss process is when three 
waves are annihilated into the vacuum. All these processes are not allowed by energy-momentum 
conservation for massive particles in fiat space. However, in dS space all these processes are allowed 



30], [3l[, [34], [331], [3j], [35|], [8|], because there is no energy conservation. 

Due to the presence of Kp in (j23p we have extra terms in CI. All of them can be obtained from 
those listed in ([25]) via the simultaneous substitutions of (1 + nfc^p)'s and nfc,p's by K^^p's or K^p's. 
E.g. we encounter terms of the form: 



[(1 + Up) Kki "-fca - np Kki (1 + rafcj] (26) 

which as well correspond to the two competing processes — one that the wave with the momentum 
p lost (gained) in such a process that instead of the creation (annihilation) of the two particles ki 
and k2 we see single k2 and the missing momentum ki is gone into (taken from) the background 
quantum state of the theory. Which should not be confused with the background geometry. 
Similarly we encounter terms of the form 



[(1 + Up) Kki - rip Kfca] (27) 

which describe the processes in which both ki and k2 are coming from (going to) the background 
state. 

Thus, we find it natural to interpret the anomalous quantum average Kp = {ap a-p) as the 
measure of the strength of the backreaction on the background quantum state of the theory of the 
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various processes described by the standard terms (p5]) in CI. Then it should be expected that in 
the vicinity of the equihbrium the backreaction is small. Which means that Kp is suppressed in 
comparison with Up and can be neglected. In such a situation the system of KE reduces to one 
equation for Up only. 

Rephrasing, if we start from an initial state, which is substantially different from the eventual 
stationary one, the generated Kp is comparable to Up. But as the state of the theory approaches 
the equilibrium, Kp becomes suppressed, which signals the small backreaction. 

The picture we have in mind is as follows. Note first of all that the in-state (say BD one) 
looks as the coherent state from the point of view of our preferable out-state (specified by the 
above out Jost harmonics). The situation is similar to the one for the QED in the constant electric 
field background [sS*]. Now consider a flat space ordinary QFT, which has the unique Poincare 
invariant vacuum state. Consider its evolution in time if the in-state is some excited coherent 
state. Intuitively it is natural to think that the final state of the theory will be build on the basis of 
the appropriate vacuum state under consideration. And one can explicitly see that in the in-state 
the anomalous quantum average {aa) will be of the same order as {a~^a), while in the out-state 
the anomalous quantum average will be suppressed. Here o"*" and a are creation and annihilation 
operators corresponding to the correct Poincare invariant vacuum state. 



IV. KINETIC EQUATION IN POINCARE PATCH 



Taking into account the one loop calculation above we conclude that if one starts from the 
BD state at past infinity then he has to solve KE equation for Up and Hp together. At the same 
time, we may assume that in the expanding PP of dS space there is a final state, which is close to 
stationarity and is build on the basis of the Fock space corresponding to the out Jost harmonics 
h(x) ~ Ji^{x). For the latter harmonics one can write the KE equation containing only rip. 

Then, as 77 — )• 0, we can put Kp = in (f23l) to arrive at 



— ^ = — / dA;r? A:r? 2 / dy' {y')^ x 
dlog{r]) Jo J 00 



X <^ Re C 



^ kii, kri, — ^ krj 
k k 



C* 



P I 1 /' 



-2Re C 



(1 + Hp) rifc - np (1 + n^) (1 + n|p_;,|) (r?) + 



, \k-p\ , P , 

kri, ki], — kri 

k k 



C* 



k y^k^ 



Uk {I + n\k_p\) {1 + Up) - {l + nk)n\k_p\np (77) + 



+Re ( D 



^ \Pj^ — j_ ^ P J 
k k 



D* 



I \p + k\ I p ; 



k y^k^ 



(l+nfc)(l+n|p+fc|)(l + np 



Uk n\p+k\ Up 



iv) } 



(28) 



where C[x, y, z] = h*{x) h{y) h(z), D[x, y, z] = h{x) h{y) h{z) and h{x) is the specified above set of 
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solutions of the Bessel equation. In deriving this equation we assumed that all quantities under 



the integral over d^k depend only on 
This is not yet the equation we are 



and changed the variables y' = krj' . 
ooking for. But before going on let us make a few comments. 
First, it is not hard to see that for short periods of time t = —\ogr] <^ 1, i.e. well within the 
cosmological horizon, and for m, /c ^ 1, the rj' integrals simplify and the above equation reduces 
to the one in flat space presented in the Appendix. However, for cosmological times one has to use 
solutions of the Bessel equation instead of the flat space plain waves. Then one does not obtain 
the ^-functions ensuring energy conservation inside the CI. That is of cause should be the case in 
such a time dependent gravitational background as dS space. 

Second, the prefactors of the Up dependent terms inside the CI can be considered as the defini- 
tions of the rates of the six processes described by the CI. This can be a way to define the rates in 
the circumstances when there is no well defined notion of the S-matrix 2C | . 



For the small physical momenta {pr] — t- 0) one should look for a solution of this equation in the 
form of the function of the physical momentum prj — np{rj) = n{prj). In fact, PP brakes some 
part of dS isometry, but is invariant at least under the simultaneous rescaling of r] and x. Hence, 
even if we start with the non-invariant under this symmetry state one can expect that due to 
the expansion, which smoothes everything out, unless backreaction becomes strong, all physical 
quantities are going to be invariant under the symmetry in question in the future infinity. In 
particular, that means that the distribution np{rf) should depend on the invariant quantity pr]. 

Indeed, for the small physical momenta the equation (|28p reduces to: 



dn{x) 

d\og{x) Air'^ fi jQ 



1 r i 
dyiy)^ / dy' {yy x 



X <^ Re <^ y 



+2Re<y 



+Re<^y 



h'' (y) 



\hiyW 



(y) 



TT e 



-TTfl 



4 sinh(7r ^u) \y\ 



TT e 



4 sinh(7r fi) \y' 
{1 + n{x)) n{y)'^ - n{x) {1 + n{y))'^ 



TT e 



-TT/l 



4 sinh(7r fi) \y\ 



/\ -in 



[y 



\Hy') 



TT e 



4 sinh(7r fi) \y'\ 



n{y) {I + n{y)) {1 + n{x)) - {I + n{y)) n{y) n{x) 



+ 



+ 



vr e 



-TTfl 



4 sinh(7r /i) \y\ 



TT e 



-TTfl 



4 sinh(7r /i) \y' 
(1 + n(i/))^ (1 + n(x)) - n{y)'^n{x) 



(29) 



where x = pr], y = kr] and y' = kij' . In the last expression we have neglected the contribution from 
the region k < p to the integral over k (i.e. over y), because the main contribution to the integrals 
comes from A; r/ ~ /i, while pr] <^ fi. 



15 



V. SOLUTION IN THE EXPANDING POINCARE PATCH 

As we explain in the Appendix in Minkowski space-time the Plankian distribution Up = 
l/(e'^/"^ — 1) annihilates the first two contributions inside the CI (I28p . The last term in the 
CI is just forbidden by the energy conservation in Minkowski space. However, neither the first two 
terms nor the last term in (j28p are annihilated by the Plankian distribution, because there is no 
any restriction on the energy in dS space. Note that even for the case Up = the CI does not 
vanish due to the particle creation out of vacuum. 

Now we assume that in the future infinity of the expanding PP, or for the low physical momenta, 
n{pr]) is very small. At the same time the particle number density for the high physical momenta 
should be even smaller, n{x) ^ n{y) for y ^ x. Hence, in this limit one can approximate 



(1 + n{x)) n{y)'^ — n{x) (1 + n{y))'^ —n{x), 
n{y) (1 + n{y)) (1 + n{x)) - (1 + n{y)) n{y) n{x) ^ n{y), 

(1 + n{y)f (1 + n{x)) - n{yf n{x) ^ 1. (30) 

Furthermore, the main contribution to the second term (I29p comes from the region y ~ ^, while 
X <C /i. Hence, we can neglect the second term on the RHS of ()29p because it is proportional to 
n{y ~ /u) ^ n{x). Then the KE reduces to: 



r 



r' 



A2 



27r2 fi 
2 7r2/i 



dyy2 



dyy" 



dn{x) 
d\og{x) 

\h''{y) ^ 



r n(x) — r', 



4 sinh(7r /i) \y\ 



IT e 



4 sinh(7r fi) \y\ 



(31) 



Which is a kind of renormalization group equation where the CI plays the role of the /3-function 
0]. In fact, we obtain an equation which shows how the distribution of particles n{pr]) changes 
with the change of the scale pij. Moreover, as follows from the discussion in the Appendix the 
solution of the KE sums up bubble diagrams, which is exactly what the renormalization group 
equations do. Similar situation is discussed in [l3 |. 
The solution of the obtained differential equation is 



r' r p ■ 



(32) 



where C is the integration constant, which may depend on the initial conditions. The obtained 
solution is valid because ^ ~ e^'^^'^ <^ 1 for /i ^ 1. Furthermore, at the leading order in we 
reproduce the one loop result of the previous section if C = —1. One should recall here that n{pr]) 
is the particle density per physical volume. 
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The obtained distribution has the stationary value e~^'^^, which approximately annihilates the 
CI in the IR limit. It is reached when the production of particles is equilibrated by their decay. 
In fact, from what we have mentioned above it is clear that T defines the decay rate of the scalar 
particle into two, while T' defines the particle production rate. Towards the future infinity logx 
is decreasing, hence, indeed V corresponds to the gain, while T to the loss in (j3ip . Furthermore, 
obviously the loss in question should be proportional to the particle density, while the gain should 
be just a constant for low distributions. 

At the stationary point in question the theory in dS space can be described by the analytical 
continuation from the sphere. But one can not describe the approach to the equilibrium ()32p via 
such an analytical continuation. 

Finally, we have obtained the solution (j32p assuming that n{pr]) decays in the future due to the 
expansion of the space and despite the constant rate of particle production. Such an assumption 
is natural if we start with low particle density at past infinity. However, there is a question if this 
situation would change or not once we will start with some different state at past infinity? Of cause 
if we will start with a very high density of particles, in comparison with the vacuum energy, then 
the situation will be much different, but what will happen say at the intermediate densities? To 
address this issue one may look for other solutions of the KE (j28p . But we find it more appropriate 
to study the problem in global dS space, because global dS sets up initial conditions for the PP — 
for the infiation. This will be done elsewhere. 



VI. SOLUTION IN CONTRACTING POINCARE PATCH 



In the contracting PP we do not expect stationary situation in the future infinity [IJ]. Hence, 
we do not expect KE to be applicable in this situation. The situation demands the consideration 
of the KE for both {a~^a) and (aa) simultaneously. But for the moment let us assume that the KE 
is applicable and see what kind of solution it leads to. We will explicitly see that for the obtained 
solution the approximation, which we have used to derive the KE, brakes down for late enough 
times. 

It is not hard to obtain the KE for the low physical momenta in the contracting patch via direct 
derivation or via the time-reversal (t -H- — t or 7/ o 1 /r]) from the equation in the expanding patch. 
The only change wrt the equation (j28p is the relative sign between the LHS and RHS. For low 
physical momenta in the contracting PP we still can use h(x) ~ Ji^{x), but now they play the role 
of the Jost functions at past infinity. 

We are interested in the behavior of the solution at the future infinity r/ — )• 00 and expect that 
for low momenta the distribution grows with time, due to the contraction of the space and constant 
particle production. As we will see in a moment the products of C[x,y,z\ and D[x,y,z\ on their 
complex conjugates in (j28p do not depend on p and k for low enough x, y and z. That is related 
to the usual flatness of the spectrum in dS background. 

Hence, it is natural to assume that np{r]) for pr] <^ fi does not depend on p and can be taken 
out from the integral over k on the RHS of (j28p . Moreover, due to the expected explosive growth 
of n{ri) for the low momenta we can make the following approximations 
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(1 + Up) rik np_k - np (1 + rife) (1 + Up^k) « -n^iv), 
Uk (1 + Uk-p) (1 + Up) - (1 + rik) Uk-p Up « n^ir]), 
(1 + nfc) (1 + Wp+fc) (1 + np) - Hk Hp+k rip « n^{r]) (33) 

if A; and p are small enough. 

The reason why the range of momenta for which our approximation is valid (prj <C /x) narrows 
down as the time goes by (r) — >■ oo) is as follows. Due to the contraction of the space in question 
the long wave length fluctuations cross into the horizon with the growth of time. But there are the 
horizon size modes which show the explosive behavior of their distribution. 

Having made these observations, let us split the integrals over y and y' in the CI into two parts: 
one due to the region where y,y' <^ n and the other due to the region where y,y' 3> Then the 
KE reduces to: 



A2 



dn{r]) 

d\og{rf) Att^ fi^ Jq 



£dy{y)^!^-Re y''" £ dy' {y')^+''' n'^{r]) + 



-2 Re 



-1 11 



+ 3Re 



Jo 



n\v) }+■■■ 



(34) 



Dots on the RHS of this equation stand for the terms which are suppressed in comparison with 
the terms explicitly written in this equation, because n{y) for y ^ /j, should be small as com- 
pared to n{ri). Thus the integro-differential KE again reduces to the differential equation of the 
renormalization group type: 



^^^^^ -Tn\r,), where T^ ^ , f,^' > 0. (35) 



d\og{r]) ' 27r2m2(m2-|) 

The solution of this equation is: 



"^^^ ^ A-tiogv ^ nkf' ^^^^ 

and is indeed independent of p. The solution under consideration is valid for rj < rjo = e'^""**/'^^ 3> 1, 
A is an integration constant, which may depend on the previous history of the evolution of the 
state and, in particular, on the initial state. 

As we see this solution has a pole at some finite r]Q. The situation is similar to the renormalization 
group Landau-Pomeranchuk pole if we consider the CI as the /3-function. In the vicinity of the 
point 770 the approximation in which we have derived the KE brakes down, but in any case we 
already have to take into account backreaction both on the state of the QFT in question and on 
the background. 
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VII. CONCLUSIONS 



Several comments are in order at this point. First, in this note we have considered 4D (f)^ QFT, 
but most of our conclusions and formulas can be easily extended to the arbitrary dimensions and 
for (f/^ theory. As well one can write the CI at higher orders in A. In fact, from the discussion 
in the main body of the text the structure and physical origin of all terms in the CI should be 
clear: one should just include into CI all possible processes, which are allowed by the momentum 
conservation, and substitute the plain waves e*'^* (inside the expressions leading to the (5-functions 
establishing energy conservation) by the out Jost harmonics. 

Second, the term in the CI (j28p . which is responsible for the creation of particles out of the 
vacuum, does not vanish as well for the odd dimensional dS space-times. Thus, there is the 
particle production as well in the odd dimensional dS space-times. This observation contradicts 
some earlier climes in the literature. 

In conclusion, we have shown that the proper interpretation of the Schwinger-Keldish one loop 
renormalization of the two point function is in terms of generation of the slow functions Up = {a^ Up) 
and Kp = (apa-p). Which signals the particle creation in dS space. We have explicitly shown that 
the one loop contribution to the two-point function does not reduce to the mass renormalization. 
Hence, loop corrections for the quantum fields in dS space can not be obtained via analytical 
continuation from the sphere. 

We have shown that for the BD harmonics the IR value of Kp is of the same order as Up. Hence, 
these harmonics are not suitable for writing KE for Up only: For BD state one has to solve the 
system of KE for both Up and Kp. We derive such a system of KE as well. However, for the Jost 
harmonics gp{r]) ~ 77'^/^ Ji/iipv) the IR value of Kp is suppressed in comparison with the one for Hp. 
We suspect that this observation means that from whatever homogeneous state at past infinity in 
PP we start eventually the state of the theory flows to the one build with the use of the future 
Jost harmonics. In fact, as we explain in the main body of the text, it is natural to interpret the 
anomalous quantum average Kp = {ap a-p) as the measure of the strength of the backreaction on 
the background quantum state of various processes described by the standard terms in CI. Thus, 
if we start from an initial state, which is substantially different from the eventual stationary one, 
the generated Kp is comparable to Up. But as the state of the theory approaches the equilibrium, 
Kp becomes suppressed, which signals the small backreaction. 

We have derived the KE containing only Up both in expanding and contracting PP of dS space 
and found solutions of these equation. These solutions can be understood on general physical 
grounds. In the expanding PP it looks like 



n{pi]) 



F' 



C {vr]f + 1 



(37) 



where C is the integration constant, which may depend on the initial conditions. The obtained 



solution is valid because 



-2 TT/i 



^ 1 for /i ^ 1. It has a stationary Gibbons-Hawking value 



e which approximately in IR limit annihilates the CI and is reached when the production of 
particles is equilibrated by their decay. In fact, from what we have stated in the main body of the 
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text it is clear that T defines the decay rate of the scalar particle into two, while T' defines the 
particle production rate. 

Note that np{r]) is the density per physical volume. The density per co-moving volume decays 



to zero in the future infinity: Up 



C [pT]) +1 (see the discussion above). 



In the contracting PP the solution for low momenta p is 



n{r]) 



1 



^ — r log rj 



f log 20 



(38) 



and is independent of p. The solution under consideration is valid for rj < r]Q = g'^""**/^^ ^> 1, yl 
is an integration constant, which may depend on the previous history of the evolution of the state 
and, in particular, on the initial state. As expected the distribution grows with time, due to the 
contraction of the space and constant particle production, and moreover has a pole at some finite 
Vo- 

The use of the expression "particle density" in this context demands some clarifications. As we 
mentioned in the main body of the text one can not diagonalize the free hamiltonian Hq in dS space 
once and forever. One can find harmonics which make Hq diagonal at past infinity, but there are 
no harmonics which diagonalize Hq in the future, because Bf^ is time dependent as r/ — )• 0. To deal 
with the appropriate notion of particles it is tempting to find such harmonics which diagonalize the 
free Hamiltonian instantaneously at a fixed moment of time r]. This can be achieved via the time 
dependent Bogolyubov transform [38|]: bk = akiij) ak + /^kir]) a'^j^, 6^ = a\{vi)al + f3l{r])a^k- 
The harmonics are 



9k{^) = algk- 1^1 gl 



dr] 



a{rj) [k"^ + a?{r]) m? 



(39) 



Hence, it is tempting to derive the KE for {b^ bp) — appropriate particle density at the given 
moment of time. 

However, because we have made a time dependent canonical (Bogolyubov) transformation to 
arrive at these harmonics we have to take into account the explicit time dependence of 6's. Then 
the problem is that the evolution equation closes wrt {bp bp) only if dtp (Sp — Up fip is negligible 
in comparison with the CI. However, this is not the case in dS space. In the latter case the 
corresponding equation does not even close wrt both (6+6) and (66). Hence, in dS space it is more 
appropriate to use the above Jost harmonics. Furthermore, if one knows {a~^a), {aa) and {a~^a~^) 
as the solutions of the KE, then he can find {b'^b), (66) and {b~^b~^) via Bogolyubov transformation 

(El. 
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IX. APPENDIX 

To make the paper self-contained in this Appendix we present the derivation of the KE for the 
flat space massless real scalar field theory with the cubic self-interaction: 

(40) 

For small A this theory describes phonons in crystals with slight non-linearity. 



S 



I 



\ (5,0)'-^ A, 



A. Operator derivation 



The general method of the derivation of the KE in the operator formalism can be found in [36l |. 
Here we concisely perform the calculations for the particular model in question. We consider the 
case when the distribution function is homogeneous in space, i.e. depends only on time, but not 
on spacial coordinates. 

The normal ordered Hamiltonian for this theory looks as: 



H = j (fkk : al at: +Hint, k = k , [uk, a+] = 6 (^k - ^ 

A f d'^kid^k2d^k3 



H.. 



int 



3 J yjki k2 ks 

3 6 [-h +k2 + h) (e-*(-'=i+'=2+^«)* : a+ a^, a^, : (-'^i+'^^+^a) * . + 

+ d(k, + k2 + ks) (e-*('=^+'=2+'=3)* : a,, a,, a^, : +e^(>^Mks)t . ^+ ^+ Al _ (4^) 



Usually one drops the last two terms in Hint due to energy-momentum conservation. However, we 
keep them to show the difference of the situation in flat space wrt the curved one. 

Consider this theory in some non-stationary initial state, characterized possibly by some density 
matrix. We would like to find how the distribution function Up = y l^a^ ap^ evolves in time. Here V 
is the volume of space, which appears because (a^ aq) = np6 {p — q) and V = 5(0). The evolution 
equation in the interacting picture is: 
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(42) 



Then: 



dup i A 



-6 y-p + ki + k2 
+2 5 (-ki + k2+P 



■ Op aki fflfca 7 e 



-i (-p+fci+fc2) t 



^ J VkTkoJ} 

i (-p+fci+fc2) i I _j_ 



g-i (-fei+fc2+p) t 



+ 5[p+ki + k2] [{-.apak.ak^ :) e 



-i (p+A;i+A;2) t 



""p ^ki ^k2 



J {-ki+k2+p) t 



J {p+ki+k2) t 



+ 



(43) 



We dropped the zero modes here. Careful study of these modes reveals the classical instability of 
the theory with cubic self-interaction. However, the terms in the CI which are responsible for this 
instability are not universal. Such terms in the cj)^ theory do not present, because one always can 
choose the vacuum value of the scalar field, which respects Z2 symmetry. 

The equation (jJS]) does not close wrt (a+ Og). So we have to find l^a'^ a+ a^) {t), (^ap a^) {t) 
and etc.. Consider e.g. 



d_ 

lit 



• ^k-^ ^k2 ^kz 



intj '■ (^ki ^kz 



(44) 



After the calculation of the commutator in this expression we see that its RHS depends on such 
quantities as e.g. ^a^^ a^^ "^fe^)- Hence, we have to derive the evolution equations for them and 
so on. This way one obtains the so called Bogolyubov hierarchy. To truncate the hierarchy one 
should decide what is the approximation he would like to consider. If one would like to find the 
CI (the RHS of (|42|) ) at the order he has to truncate the sequence of these equations already in 
). This can be done via the application of the Wick's theorem, e.g.: 



Ofci <4^ flfca a^) = (flfci flfe^) (ofc2 afca) + (ofci a^^) (afcg ^4) 



(45) 



and so on. The justification of such an anzats/approximation is given in the next subsection of the 
Appendix (see as well |37l|). 



So, if one uses the Wick's theorem in the commutator 



Hintt ^k's 



he obtains: 



d 

^ V "fei "^=2 '^^^ ■ 



5 (-ki + k2 + h) e^(-^i+'=2+'=3)* 



i2A 

(1 + rifcj rifcj n^^ - n^^ (1 + n^-J (1 + n^g) 



where apart from (a^ a^) = np6 {p — q) we have used that (a^ a+) = (1 + Up) 6 {p 
we can find the equation for ( a^^ a^^ ^fca / • result: 



(46) 

). Similarly 
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-i ( — fcl+fc2+fc3) t 



(-fcl+fc2+fc3) t 



■ ofei % % 7 W 



: Ofcj -.j (t) - e 

^ 2 ^ S (-k, + ^2 + ^3) Re r dt'e-* (-fci+fe+fo) {t-t') ^ 



it') 



(1 +nk^)nk^ rij^^ - (1 + n^J (1 + 77.^3) 



(47) 



Here to is the moment of time when we switch on interactions. 

Along the same hnes, one can find the answer for e"*^'^^^*^^^'^^^* (: Ofci Ofe2 t^fcs '■ 
gi(fei+fc2+fc3)t /. ^+ ^+ ^+ .\ 'pj^Tjg^ -^g obtain the following equation: 



P 



drip 
It 



ki k2 



(1 + np) nfci - np (1 + ) (1 + 



S \^-p + ki + A;2j J dt' cos [(-p + fei + ^2) (t - t')] 

+2 5 (^-ki + k2 + J^^ dt' cos [{-ki + k2 + p) (t - t')] n^^ (1 + n^^) (1 + rip) - (1 + J np (t') 

+ S (jci + k2+pj dt' cos [{ki + k2+p) (t-t')] (1 + n^J (1 + n^J (1 + Up) -nk^Uk^Up (t')|(48) 

We assume that n^'s depend on t slowly and can be taken out from the integral over t' . Then we 
take t — )• +00 and to ~^ —00. The result of the integration over t' is the minus (5-function ensuring 
the energy conservation in the corresponding process. Hence, the last term in (j48p vanishes because 
it describes processes forbidden by the energy-momentum conservation. In the case of massless 
theory the first two terms describe allowed processes and the KE acquires the form: 



P 



drir, , , 9 f d^ki d^k2 
— 4 A / — ^ — ^ X 



dt 



ki k2 



5 ( -p+ A;i + /c2 I 5 {-p + ki + h 



(1 + Up) Uk^ rif^^ - np (1 + rifcj (1 + n^J 

+ 25(^-ki + k2+^ 5{-ki + k2+p) n^.^ (1 + n^J (1 + np) - (1 + n^J np | . (49) 

Having understood the above equation and the physical meaning of all terms in it (see the discussion 
in the main body of the text) it is not hard to write the KE for the massive real scalar (j)'^ theory: 
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dnp 
~dt 



UJi UJ2 



X 



3 5 ( -p- ki+k2 + k^ 



sin [{-uj - UJ1+U2 + UJ3) {t - to)] 



-UJ - UJi + UJ2 + UJ3 

(1 + n^^) (1 + 

(1 + 



X 



+3(5 -fci + A;2 + /i;3 +p 



sin [{-uji + UJ2 + UJ3 + (jj) {t - to)] 



null (1 



+6 



-uJi + a;2 + W3 + 
n^^J (1 + n^^3)(l + n^^) -(l + n^Jn^, 

sin [(-tj + wi + a;2 + ws) {t - to)] 



X 



-p + /Cl + ^2 + ^3 

+ 



+(5 /i:i + A;2 + /i:3+P 



-w + wi + a;2 + t^s 

nui (1 + (1 + »1-a;2)(l + "-^a) 

\ sin [{u!i + u!2 + UJ3 + uj) {t - to)] 



OJi + UJ2 + U}3 + ^ 



(1 + Uuii) (1 + n^^a) (1 + "-t^s) (1 + "-a;) 



(50) 



00 the function is reduced to the 5- 



where = + m^. As t — to 
functions ensuring energy conservation. The only allowed by the energy-momentum con- 
servation process is the scattering between the scalar particles, i.e. the CI contains only 
6 {ui + UJ2 - - ^) [nuji nuj2 (1 + ^^^3) (1 + n^j) - (1 + n^^i) (1 + n^^a) "-^3 n-uj]- Two more terms 
will stay in the CI for the massless cj)'^ theory: only the last term in the above CI is forbidden 
by the energy-momentum conservation in the massless theory. 

In the small density limit, n^j <^ 1, for all u we can neglect n^j in comparison with 1 in all 
expressions inside the CI^. Then for the massive c/)^ theory the CI will reduce to its appropriate 



classical Boltzmann form 6 {coi + u;2 — ^3 — a;) [n^^-^ 



Now one can immediately see 



that because of the energy conservation this expression vanishes for the equilibrium Boltzmann 
distribution n^j oc e~T ^ for some constant T. As well it is not hard to see that all terms in the CI 
(HH) and (f50|) (for {t — to) ^ 00) vanish when n^^ = ^ui/t_i ■ 

Finally, it is not hard to generalize the KE to the spatially inhomogeneous case by the calculation 
of d^n = y {[Pf^,a~^ a]), where = {H,Pi) is the momentum operator. Then the LHS of (|49p 
and (|50p will change to d^np{x). The RHS will be the same. 



B. Derivation from the Dyson-Schwinger equation 



The general method to derive the KE from the 
Keldysh) formalism can be found in e.g. 2j] or 



^ 3SE of the in-in (non-stationary or Schwinger- 
25( 1 . Here we concisely repeat that derivation. 



* This is, in fact, proper classical limit, because otherwise due to high particle density quantum coherence between 
particles is not destroyed. 
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The main reason to present this derivation here is to show the relation of the KE to the summa- 
tion of the loop diagrams in the IR limit. The technical explanation of the relation between the 
Wick contractions, which were used in the derivation of the Boltzmann equation, and the partial 
resummation of loop graphs can be found in [37( . 

The reason for the Schwinger-Keldysh diagrammatic technic can be seen from the following 
observation. Suppose we would like to calculate the expectation value of an operator O at some 
moment of time t wrt some state ^ 



(O) it) 



(51) 



where H{t) is the full Hamiltonian of the system, which may depend on time due to the presence of 
e.g. time dependent background fields, T is anti-time ordering, to is initial moment of time when 
we know the expectation value. Transferring to the interaction picture, we get: 



{O) it) = {^\S+it,to)T[Oo{t)Sit,to)]\^) = 
= \S+{t,to)S+{+oo,t) S{+oo,t)T poit) S{t,to)]\^) = 

= {^\S+{+oo,to)T[Oo{t)S{+oo,to)]\^), (52) 

where S{t,tQ) = Te * /*(, '^^ ) and C'o(t) is the time dependence of the operator O in the 
interaction picture. Adiabatic switching off of the interactions is assumed in the future infinity. 

If we adiabatically switch on the interaction around the time Iq, then we can write the expec- 
tation value as: 



(O) (t) = |5+(+oo, -oo) T [Oo{t) S{+oo, -oo)]| ^> . (53) 

Good question is that if one can take to to — oo. If the state of the theory does become stationary 
(e.g. thermalizes), then to can be taken to the past infinity. However, if the state does not get 
stationary, which may be the case if the background field is never switched off, then to can not 
be taken to minus infinity because of the explosive behavior of the correlation functions when 

(t — to) —5- OO. 

Now if |\I') were the true vacuum state \vacuum) of the free theory, then by adi- 
abatic switching on and off the interactions one could not disturb such a state, i.e. 
{vacuum |5^(+oo, — oo)| excited state) = 0, while {{vacuum |5'^(+oo, — oo)| vacuum)\ = 1 and 



(O) (t) = (vacuum 1 5^ (+00, —oo) I state) (state |r [Oo(t) 5(+oo, —oo)]| t;acuum) = 

state 

= (^vacuum |S^(+oo, — oo)| vacuum) {vacuum \T [C'o(t) 5(+oo, — oo)]| vacuum) = 

{vacuum \T [C'o(t) S'(+oo, — oo)]| vacuum) 
{vacuum |S'(+oo, — oo)| vacuum) 
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This way we arrive at the standard Feynman diagrammatic technic, which is based only on the 
T-ordered quantities. 

However, if l^') is not a stable state one can not use the above machinery and has to deal directly 
with ()53p . One has to expand both S and in powers of the interaction Hamiltonian and apply 
the Wick's theorem. Then one will encounter two types of vertices and four types of propagators. 
The vertices will be coming from S and . At the same time propagators appear from the Wick 
contractions inside S (time ordered), from those inside (anti-time ordered) and from Wick 
contractions between S and S^. However, there are only three independent propagators. 

The functional integral in such a situation has the action, which schematically can be written 



as 



25|: 



S= I dt I d^x 



(55) 



where C is the Keldysh time contour running from — oo to +00 and back. That is due to the 
presence of both S and S'^ in ()53p . The exact expression for the action will be given in a moment. 
This action can be rewritten as: 



S 



+00 



dt / d^x 



(56) 



where is the field on the direct part of the time contour, while 0_ is the field on the reverse 
part of it. The kinetic term in this equation again is written schematically 25|. After the Keldysh 
rotation of the fields: 



V2 



the precise form of the action is as follows 25l |: 



V2 



(57) 



S 



d^xfdSiM^) ,Mx))(^ [Dq-\^.y)\fMy) 

J yyy^ciy )j y [D^]-\x,y) [D,Y i^^v) J \ Mv) 

1 



2 A d^x 



(58) 



The vertices and propagators in this theory are shown in the figure. Here Dq and Dq are the 
advanced and retarded propagators, whose Fourier transforms look as 



{e±iOy - fc2 



(59) 



In the x-space they have the form 
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Figure 3: The solid line corresponds to (/)c;, while the dashed line — to (f)q. 



= 0(±At)5(^)(x-y)n, (60) 

where □ is the D'Alamber's operator. The retarded and advanced propagators carry information 
about the spectrum of the theory and have the following relevant for us properties: 



D^ix,y)D§ix,y) = 0, D^it,t) + D^it,t) = 0, [D^f = D^, (61) 

which remain valid even for their quantum corrected versions ^2^. 

The Keldysh propagator [^o^]"*" = ~Dq carries statistical information about the the- 
ory: it shows which levels from the spectrum are occupied. By definition D(f{LO,k) = 
/ d'^xe^^'^^^ ({(/)(x), (/)(0)}), where {•,•} means the anti-commutator. For the thermal state it 
acquires the form: 



2T 



coth^ [D§{uj,k) - D^{uj,k)] = coth ^ 5 



2T 



UJ 



k' 



(62) 



This propagator is present in the above action to regularize the functional integral 25l |. 

The last expression allows to guess the ansatz for the full quantum Keldysh propagator 25|, 
when one is close to the stationary situation: 



Z^^(x, v) = j d^z [D^{x, z) F{z, v) - F{x, z) D^{z, y)] ^ [D^ o F - F o D^] {x, y) (63) 

with some unknown in general kernel F{x, y) which characterizes the statistical distribution. In 
the stationary situation the Fourier transform of F{x,y) is coth ^ = 1 + 2n^, where n^^ is the 
Plankian distribution function. Furthermore, 
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[^0 



11^ 



[D^]-'oD^o[D^]-'=[D^ 



J 



(64) 



One can write the DSE for the full quantum matrix propagator: 



D 







which, as can be shown 



251], keeps the same form as Dq. The equation looks as: 



(65) 



D = h 



where 



(66) 







(67) 



is the self-energy matrix. It can be shown that it has the same form and properties as Dq ^ . Below 
we will see this fact at one loop level. In components of D and S the DSE can be written as 



[□,F] = S^- (S^oF-FoS^) (68) 

if one uses the ansatz = D^o F — Fo and neglects [-Cq hi comparison with S^, because 
it is just a regulator. 

Let us derive explicitly the DSE in one loop approximation. Note that while considering full 
we take the bare D^'^. The justification of this approximation comes from the observation 
that the only way renormalization appears in and is through the change of the spectrum 
— renormalization of the mass and etc.. So once we know the spectrum of quasiparticles precisely 
it means that we know the retarded and advanced propagators as classical objects. As well we 
substitute into the DSE the IR, i.e. renormalized, value of the vertex in the theory. In such cir- 
cumstances the DSE becomes an equation for only — for the propagator containing statistical 
information about the state in the theory. 

First we calculate the self-energy at one loop order: 

• It is straightforward to see that S^-d is zero (see fig.): 

T.,i_,i^\^D^{x,y)D^{x,y)=0, (69) 
because of the presented above properties of and Dq. 
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Figure 4: The graph defining T,cl-cl- 




Figure 5: The graph defining S"^. 

• At the same time (see fig.) 

Ee/-g = S^(x,y) =4iA2l)^(x,y)D^(x,y) /O (70) 

Since T,^ ~ ~ 9{ty — tx), this quantity is indeed of an advanced type and should, as 
it does, stand in the upper triangular corner of the S matrix. The prefactor 4 comes from 
4 ways of choosing external legs, 2 internal permutations and 1/2! for having two identical 
vertexes. 
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Figure 6: The graph defining Y,^ . 

Similarly (see fig.) 

= S^(rc, y) = AiX^ D§{x, y) D^'ix, y) (71) 
As well ~ ~ e{t^ - ty). Hence, at the one loop level = [S^]^. 
And finally (see fig.): 



^ E^(x, 2/) = 2 [D^{x, y)f + Gi^\[D\x,yf+Qi\^ [D\x, y)] ' 



2iA2 ([Z)^(x,y)]' + [D^(x,y) -D^(a;,2/)]') , (72) 



where at the last step we have used the property D^{x,y) D^(x,y) = 0. Now one can see 
that because [D^]^ = and [D^]^ = -D^ we have that [S^]'^ = -S^. 

Thus, as promised at one loop level S has the same properties as Dq^- Now we plague these 
expressions for the components of S into the DSE for F and use the Wigner transformation: 



A{x,y) = A{X,x) , X = ^^—^, x = x-y 

A {X, x) = j (fke' ''"'^^ a {X, k) . (73) 

Here a is the Wigner transform of A. 



The Wigner transformation has the following properties 25|]. If 
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A {X, x) = j d^k e' ^'^^^ a {X, k) , B{X,x) = j d^k e' ^'^^^ b {X, k) , (74) 
and if C{x, y) = A{x, y) B{x, y), then 

c{X,k) = j d\e-'^''^>'C{X,x) = J d%d'^k2S^^\k-ki-k2)a{X,ki)b{X,k2). (75) 
At the same time if C = Ao B, then: 



:{X,k)=a{X,k) e '(^^ b{X,k) ^ a{X,k) b{X,k) + ... 



(76) 



The last approximation is vahd if a and b are relatively slow functions of X = {x + y)/2 and 
fast functions of % = a; — y. Thus, we should have a separation of scales in the problem under 
consideration. 

Using this approximation, the derived above expressions for the self-energy and defining as d^, 
d^ and / the Wigner transforms of Dq, Dq and F, correspondingly, we obtain that: 

T.^ o F - F oT,^ ^ f {X,k) j d'^ki d^fcs S^'^'^ {k-ki- fcs) {X, ki) d^ {X, k2) ~ 

2 i AV (X, k) j d^ki d''k2 (J(^) (A; - /ci - k2) Ad {X, ki) Aa {X, ^2) [/ {X, /C2) + / {X, ki)] (77) 

where = [d^ — d'^] and on the last step we have substituted the expression for d^ through d^, 
d^ and /. At the end we performed the symmetrization of the argument of / under exchange of 
ki <H- k2. 
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Similarly: 



^ 2 i A2 y" d^ki d^k2 <5(4) {k-ki- k2) Aa {X; ki) (X; fca) [/ (X; ki) f (X; ^2) + 1] . (78) 



And finally, 



[F,a]^-2ik^^f{X;k) 



Putting all this together, we obtain: 



(79) 



F df, f {X, k) = 2X'^ j (fki (fk2 5^^) {k-ki- k2) Ad{X, ki) Ad{X, yts) x 

x[f{X,ki) fiX,k2) +l-fiX,k)[fiX,ki) + f{X,k2)]} (80) 
Now due to the properties of the propagator , the Wigner transform of F obeys: 



fiX,k) = -fiX,-k) (81) 

Furthermore because A^(A;) ~ DQ(k) — DQ{k) ~ 6 (Jvq — k^^ the ki and k2 legs are on mass-shell. 
We as well put k on mass-shell — ko = ± k . Then the mass-shell distribution function obeys: 



f {X,k) = signiko) f ( X, sign{ko) k 



(82) 



Representing the / function through the distribution function, f {X;k) = 1 + 2n{X;k), we see 
that equation ([80]) reduces to the spatially inhomogeneous form of the KE obtained in the previous 
section of the Appendix. 



Note that the whole CI comes from imaginary contribution to S |25l ]. Apart from that, the 
positive contribution to the CI (the gain processes) comes from S^, while the negative (the loss) 
one comes from o F — F o S"^. 



[1] E. Mottola, Phys. Rev. D 31, 754 (1985). 

I. Antoniadis, P. O. Mazur and E. Mottola, New J. Phys. 9, 11 (2007) ^arXiv:gr-qc/0612068| ; 

E. Mottola, Phys. Rev. D 33, 1616 (1986); 

E. Mottola, Phys. Rev. D 33, 2136 (1986); 

P. Mazur and E. Mottola, Nucl. Phys. B 278, 694 (1986); 

P. R. Anderson, C. Molina-Paris, E. Mottola, Phys. Rev. D80, 084005 (2009). |arXiv:0907.0823l [gr-gc]] : 
E. Mottola, farXiv: 1008. '50061 [gr-gc]] . 
[2] B. Allen, "Vacuum States In De Sitter Space," Phys. Rev. D 32, 3136 (1985). 



32 



[3] A. D. Dolgov, M. B. Einhorn and V. I. Zakharov, Phys. Rev. D 52, 717 (1995) [arXiv:gr-qc/9403056] ; 

Acta Phys. Polon. B 26, 65 (1995) arXiv^r-qc/9405026'. 
[4] N. C. Tsamis and R. P. Woodard, Class. Quant. Grav. 11, 2969 (1994); 

N. C. Tsamis and R. P. Woodard, Phys. Lett. B 292, 269 (1992); 

N. C. Tsamis and R. P. Woodard, Commun. Math. Phys. 162, 217 (1994); 

N. C. Tsamis and R. P. Woodard, Annals Phys. 238, 1 (1995); 

N. C. Tsamis and R. P. Woodard, Nucl. Phys. B 474, 235 (1996) arXiv:hep-ph/9602315"; 

N. C. Tsamis and R. P. Woodard, Phys. Rev. D 54, 2621 (1996) arXiv:hep-ph/96023l7 ; 

N. C. Tsamis and R. P. Woodard, Annals Phys. 253, 1 (1997) | arXiv:hep- ph/96023l6 . 
[5] J. Garriga and T. Tanaka, Phys. Rev. D 77, 024021 (2008) [arXiv:hep-th/0706.0295]. 
[6] S. Weinberg, Phys. Rev. D 72, 043514 (2005) |arXiv:hep-th/0506236| ; 

S. Weinberg, Phys. Rev. D 74, 023508 (2006) |arXiv: hep-th/0605244 . 
[7] A. M. Polyakov, Nucl. Phys. B 797, 199 (2008) [a?Xiv:0709.2899 [hep-th]]. 
[8] A. M. Polyakov, Nucl. Phys. B 834, 316 (2010) [arXiv:0912.5503 [hep-th]]. 

[9] G. Perez-Nadal, A. Roura and E. Verdaguer, Class. Quant. Grav. 25, 154013 (2008) [arXiv:gr- 
qc/0806.2634]; 

Phys. Rev. D 77, 124033 (2008) [arXiv:0712.2282 [gr-qc]]; 

G. Perez-Nadal, A. Roura and E. Verdaguer, JCAP 1005, 036 (2010) [ arXiv:0911.4870l [gr-qc]]: 

A. Roura and E. Verdaguer, Int. J. Theor. Phys. 38, 3123 (1999) [arXiv:gr-qc/9904039] . 
[10] A. Higuchi, Class. Quant. Grav. 26, 072001 (2009) |arXiv:0809. 12551 [gr-ac]]: 

A. Higuchi and Y. C. Lee, Class. Quant. Grav. 26, 135019 (2009) |arXiv:0903.388Tl [gr-qc]]. 
[11] E. Alvarez and R. Vidal, JHEP 0910, 045 (2009) ^arXiv:0907.2375l [hep-th]]: 

E. Alvarez and R. Vidal, .arXiv: 1004.4867 [hep-th]. 
[12] C. P. Burgess, R. Holman, L. Leblond and S. Shandera, larXiv:1005.3551l [hep-th]. 

C. P. Burgess, L. Leblond, R. Holman and S. Shandera, JCAP 1003, 033 (2010) |arXiv:0912.T608l 
[hep-th]]. 

[13] D. Marolf and I. A. Morrison, [ar2CLv: 1006.0035' [gr-qc]; 

D. Marolf, I. A. Morrison, grXmIl04.4343 [gr-qc]]; 

A. Higuchi, D. Marolf, I. A. Morrison, Phys. Rev. D83, 084029 (2011). |arXiv: 1012.34151 [gr-qc]]: 

D. Marolf, I. A. Morrison, [ arXiv:1010!5327l [gr-qc]]. 

[14] D. Krotov, A. M. Polyakov, Nucl. Phys. B849, 410-432 (2011). |arXiv: 1012. 21071 [hep-th]]. 

[15] P. .Burda, JETP Lett. 93, 632-637 (2011). arXiv:110L2624| [hep-th]]. 

[16] S. B. Giddings, M. S. Sloth, JCAP 1007, 015 (2010). larXivil005.3287 [hep-th]]; 

S. B. Giddings, M. S. Sloth, JCAP 1101, 023 (2011). |arXiv:10 05.1056 [hep-th]]; 

A. Riotto, M. S. Sloth, JCAP 0804, 030 (2008). larXiv:0801. 18451 [hep-ph]]. 
[17] V. K. Onemh, R. P. Woodard, Class. Quant. Grav. 19, 4607 (2002). |gr-qc/020406^ ; 

V. K. Onemli, R. P. Woodard, Phys. Rev. D70, 107301 (2004). gr-qc/0406098'; 

E. O. Kahya, V. K. Onemh, Phys. Rev. D76, 043512 (2007). [gr-qc/0612026,; 

E. O. Kahya, V. K. Onemh, R. P. Woodard, Phys. Rev. D81, 023508 (2010). |arXiv:0904.48TTl [gr-ac]]. 
[18] W. Xue, K. Dasgupta, R. Brandenberger, Phys. Rev. D83, 083520 (2011). [ arXiv:1103.0285l [hep-th]]. 
[19] A. Mironov, A. Morozov, T. N. Tomaras, [arXiv:1108.2821 [gr-qc]]. 

[20] E. T. Akhmedov and P. V. Buividovich, "Interacting Field Theories in de Sitter Space are Non-Unitary," 
Phys. Rev. D 78, 104005 (2008) [ar Xiy_L0_ 808,4106 [hep-th]]; 

E. T. Akhmedov and P. Burda, "A simple way to take into account back reaction on pair creation," 
larXiv:0912.3435l [hep-th] ; 

E. T. Akhmedov, "Real or Imaginary? (On pair creation in de Sitter space)," larXiv:0909.3722l [hep-th] : 



33 



E. T. Akhmedov and E. T. Musaev, "Comments on QED with background electric fields," New J. 

Phys. 11, 103048 (2009) arXiv:090TI)424 [hep-ph]]. 
[21] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B 724, 295 (2005) arXiv:gr-qc/05051]J. 
[22] A. A. Starobinsky, J. Yokoyama, Phys. Rev. D50, 6357-6368 (1994). astro-ph/940701^ . 
[23] M. van der Meulen, J. Smit, JCAP 0711, 023 (2007). |arXiv:0707.0842l [hep-th]]. 
[24] L. D. Landau and E. M. Lifshitz, Vol. 10 (Pergamon Press, Oxford, 1975). 

[25] A.Kamenev, Book to appear. The Book corrects the collision integral for the phonons derived in the 

following review, "Many-body theory of non-equilibrium systems", arXiv:cond-m at / 041 2296^ 
[26] E. Calzctta, B. L. Hu, Phys. Rev. D35, 495 (1987); 

E. Calzctta, B. L. Hu, Phys. Rev. D37, 2878 (1988); 

E. Calzctta, S. Habib, B. L. Hu, Phys. Rev. D37, 2901 (1988); 

E. Calzctta, B. L. Hu, Phys. Rev. D40, 656-659 (1989); 

E. Calzctta, B. L. Hu, Phys. Rev. D40, 380-389 (1989). 
[27] H. Kitamoto, Y. Kitazawa, Nucl. Phys. B839, 552-579 (2010). !arXiv:1004.2451l [hep-th]]. 
[28] T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. Lond. A 360, 117 (1978). 
[29] D. P. Jatkar, L. Leblond, A. Rajaraman, [arXiv:ll07.35r3l [hep-th]] . 
[30] N. Myrhvold, Phys. Rev. D 28 2439 (1983). 

[31] D. Boyanovsky, H. J. de Vega and N. G. Sanchez, Phys. Rev. D 71, 023509 (2005) 
| arXiv:astro-ph/04094 06 ; 

D. Boyanovsky and H. J. de Vega, Phys. Rev. D 70, 063508 (2004) | arXiv:astro-ph/04 0628';^. 
[32] J. Bros, H. Epstein, M. Gaudin, U. Moschella and V. Pasquier, "Triangular invariants, three-point 
functions and particle stability on the de Sitter universe," Commun. Math. Phys. 295, 261 (2010) 
[arXiv:0901.4223 [hep-th]]; J. Bros, H. Epstein and U. Moschella, "Particle decays and stability on the 
de Sitter universe," larXiv:0812.3513l [hep-th]; 

J. Bros, H. Epstein and U. Moschella, "Lifetime of a massive particle in a de Sitter universe," JCAP 

0802, 003 (2008) | arXiv:he p-th/0612184 . 
[33] G. E. Volovik, [arXiv:gr-qc/0803.3367]; JETP Lett. 90, 1 (2009) [arXiv:gr-qc/0905.4639]. 
[34] E. T. Akhmedov, P. V. Buividovich and D. A. Singleton, arXiv:0905.2742 [gr-qc]. 
[35] E. T. Akhmedov, A. Roura and A. Sadofyev, Phys. Rev. D 82, 044035 (2010) [ar Xiv:1006.32 7l [gr-qc]]. 
[36] D.N.Zubarev, V.G. Morozov and G.Ropke, "Statistical mechanics of non-equilibrium processes", 

Akademie Verlag (Berlin), 1996. 
[37] R.D.Mattuck, Appendix of "A guide to Feynman diagrams in the many-body problem" , 2nd Edition 

(McGraw Hih, Inc., 1976). 

[38] A. A. Crib, S.G. Mamaev, V.M. Mostepanenko, Quantum effects in strong external fields, Moscow, 
(1980). 



